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Нелінійна гранична задача Алексідзе для рівняння Лапласа постулює аналітичне продовження сили тя- 
жіння в глобальних областях. Вказано алгоритм її розв'язання - визначення густини простого шару з інте- 
грального рівняння Фредгольма з регуляризацією підінтегральних похідних збурювального потенціалу. 


Тре попіїпеаг боцпаагу АЇехідге ргобіет їог І аріасе'я едиайоп розішіагез Фе єгамігу апаїунса! ргоїопрайі- 
оп іп 8ЇоБа! агеаз8. Ап аїдогійт ої 1ї5 50Їийоп 15 роїпіеа оиї аз а саїсціайоп ої 5ітріе Іауег Феп5ігу їгога ре Ртгед- 
роїт іпіестаї едцайоп мій Бе гериіагігацйоп ої зибіпіеєта! дегіуайуєез ої Ше дїзішгріпе рогепиа!. 


Нагальність вирішення нелінійної граничної задачі Алексідзе постульована непридатністю на- 
явних схем відновлення потенціалу в глобальній області. В рамках геогустинного моделювання 
регіональних структур слід розв'язувати задачу аналітичного продовження аномалій сили тяжіння 
| Дубовенко, 20094). 

Нелінійна гранична задача Алексідзе | Дубовенко, 2009а; Дубовенко, 20038| для рівняння Лап- 


ласа означає - знайти функцію У), х є у", яка задовільняє всередині необмеженої замкнутої 


області у - у (ду рівнянню Лапласа ЛУ (ж) є 0, хє у", ав точках ляпуновської межі бу облас- 





2 
ті 1ї на нескінченності - умовам: менай - ах), х є бу, У(х)-»0, х| -200, де «(х) з - МУ(о) -з 


СО дх 








з (2 (х), п (х)) - задана неперервна функція, у - обмежена область з масами Зем- 


лі, у" - її необмежене доповнення без мас, ду- межа у |і у" на поверхні Землі. 

Гармонічну в області у" функцію Ух), х є у" відшукують як потенціал простого шару |Чор- 
ний, 1995; Дубовенко, 2010| з невідомою густиною с), х є ду (інтегрованою, а загалом більш 
гладкою), поширеною на поверхні Ляпунова ду. Рівняння, з якого відновлюють невідому густину 
сх) за заданими на поверхні ду значеннями модуля градієнта потенціалу сили тяжіння (МГПСТ) 
«(х) , виведено, виходячи із аналітичного зображення сили тяжіння | ДДубовенко, 20096): 

















з"(х)е а ) а 1 со5(р, 445, 5, й б, 
де одиничні вектори р 1 4, спрямовані відповідно з точки х в точки 2 1 7), що пробігають повер- 
хо -с. х -Т; з 
хнею бу, мають вигляд р; з сов(р,х,)- ; 2 , да со5(х, 4) - ні а кут між векторами -- 
соз(р,д)- Учаах авісд) а У оз Хо 
гаю й ічІ оо оо 








Аналітичні властивості функції МГПСТ «(х) характеризуються негармонічністю | Алексидзе, 


19851, неперервністю на поверхні Ляпунова |ДДубовенко, 2009а| 1 розривністю при її перетині, про 
що свідчать такі твердження. 


Лема 1. МГПСТ не задовільняє рівнянню Лапласа в жодній точці області у . 

Лема 2. МГПСТ простого шару неперервний для будь-якої точки х є ду, яка рухається повер- 
хнею бу Ляпунова. 

Лема 3. МГПСТ простого шару, поширеного на сфері радіуса р з одиничною поверхневою 
густиною с(х) - 1, хє ду, дорівнює 
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Отже, функція а") розривна, має розрив неперервності при переході точки х через поверх- 
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ек р) дз. 
47 5, | хо 4 





ню ду. Величина цього розриву при с(х)єії, хєдбу за умови а и для сфери 


дорівнює 2. (х)- 8: (х)- 3/4, х є ду ГДубовенко, 20096). 
Граничні дані задачі відновлення потенціалу за значеннями МГПСТ «(х) а Е- МУ.) у рамках 


прийнятої моделі | Дубовенко, 20096) описує вираз 
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За припущень |Дубовенко, 2009а| про нормальний потенціал (/ (х) і потенціал притягання 
Ух) 1 (х)- Т(х) для розв'язання задачі Алексідзе запропоновано такий узагальнений ітерацій- 
ний алгоритм: 

1. за знайденими попередньо наближеннями со8(п,Х, уз 2 ом), хеду, К -« 1,2,3 напрям- 


них косинусів нормалі п (х) визначаємо на межі бу за формулою (2) і -- 1-ше наближення сили тя- 


клани У зави)- зго д9 7| РР 


Кеї дх ї 


жіння 


3 
сок(п,, т)- У сов(п, ух, )сов(ху, т), хеду, ФС )- ня (х)сов(п,, т) - у(х)сов(у, т) хедбу; 
Ка 
2. знаходимо розв'язок зовнішньої задачі Неймана для рівняння Лапласа 


мачо, хе» СТ) а) кед СО 





-»ю, (3) 


як потенціал простого шару мас "бо с б У, х єду, розподілених на поверхні ду: 


Тл Се Левине 4) 
3. невідому густину я з нелінійного інтегрального рівняння Фредгольма 2-го роду 
ба (х уч | к Х, 2Б (С м5. - -2Ф,а(х ) хе бу з () 
РЕП о ДР 1 М т) - 


де К(х,є)-е 








: соки, т) -Уобльх й ,ивхоб; при цьому 


2л дт, х-А 2 ж - 4 -б 
виходимо з постулату дея 1976), що потенціал лише тоді и розв'язком лінійного еліптич- 
ного рівняння 2-го порядку, коли його густина підкоряється інтегральному рівнянню Фредгольма 








2-го роду, яке завжди має єдиний розв'язок в області у досить малої міри. Параметри цієї області 
визначають розв'язність задачі Алексідзе, а, відтак, і роздільну здатність (компетенцію) засновано- 
го на ній методу, тому її визначення заслуговує окремої статті. 

4. Розв'язавши рівняння (35), наближено обчислимо з використанням (4) похідні 


(іча) 
уд) я ОУИ (х) 2 20(х й ЛЕ х) ше х). 1 Хі За, 5, С. Ре зай 25 (6) 
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після чого знову повторюємо алгоритм, починаючи з 1-го пункту. 


Зазначено ГДубовенко, 2009а|, що похідні (6) визначаються у внутрішніх точках області у ,а 
на межі ду їх неможливо знайти через неінтегровні особливості у підінтегральних функціях. Для 
здійснення обчислень слід знати значення похідних збурювального потенціалу саме на межі ду,і 
для їх обчислення потрібна спеціальна регуляризація інтегралів (6). 

Загалом задача обчислення вищих похідних потенціальних функцій є некоректною і має два 
основні напрямки |Чернь/й, Гольцев, 1980): 

1. виведення інтегральних зображень гармонічних функцій за значеннями її похідних на межі, 
що зводиться до вияснення умов диференційовності інтегралів, залежних від параметра (коорди- 
нат профілю) та їх чисельного інтегрування 

2. вивчення диференціальних властивостей гармонічних функцій методами кінцевих різниць. 

Регуляризація розбіжних інтегралів типу (6) зводиться до обчисленя коефіцієнтів у розкладан- 
ні за сферичними функціями за допомогою формули Гріна з граничним переходом від обмеженої 
замкнутої області у до необмеженої області у , де неперервність підінтегральної функції пору- 
шено. При цьому спочатку знаходять ті члени розкладання, які зумовлюють розбіжність поверхне- 
вих інтегралів, а потім їх віднімають з підінтегральної функції, після чого обчислюють остаточні 
значення похідних у заданій точці. При цьому отримані невласні інтеграли існують лише у смислі 
головного значення. 

ОТ, (х) 


В застосуванні до (6) вищесказане означає, що потрібно функцію Бена замінити різницею 
Хо. 
/ 
ОТ а (2,0) ОТ (2.0) 
Ох, дх; 
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, яка забезпечить збіжність інтегралу при 2 - х 1 не порушить її на нескін- 


ОТ (2) 
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Ко 


ченності (причому перший доданок задається відрізком розкладання в ряд Тейлора. Усі 
способи регуляризації розбіжних інтегралів (фактично - побудови деяких узагальнених функцій) 
відрізняються лише сингулярними функціоналами, зосередженими в особливих точках інтегрованої 
функції. Щоб уникнути розбіжності інтегралів типу (6), їх розбивають на дві частини, причому ін- 
тегрування на за інтегралі здійснюють за правилом граничного переходу: 


с, І - йт має їж хх, - 2 йт ще йт | є га 2 - 2 ї і 
Ох, | Ух ем (2-х) но о) М 
Крім того, для прискорення збіжності інтегралів можна спеціально підібрати ядро оператора 
рівняння (6) на основі компактного розщеплення |Черньг6й, Якимчик, 2001), яке є, по-суті, сумою 
двох ядер (ядра скінченного рангу та малого ядра); це розщеплення застосовне до рівнянь з ком- 
пактним оператором. 











ча че 
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Заміна коректної задачі (3) розв'язком граничної задачі (4-5) і збіжності наближень УЗ) до 
потенціалу притягання Ух), хє у" обгрунтувана в |Дубовенко, 2009а, Дубовенко, 2010). За не- 
перервної на межі ду функції густини потенціалу простого шару (6) граничні значення частинних 
похідних потенціалу 1-го рн б | Дубовенко, 2010| 








от (х Х 1 
уні ее і б їі (2М8у 5 бно м, 
але для практичного реве зрмнни варто використовувати еквівалентну формулу 
от(х х,- 1 
и Пе (од) «реся 48 5 бо) совбху ть). 


Це нелінійне рівняння за сій геометрії контактної поверхні ду(у) на класі Ляпунова 


С у) стає лінійним і має однозначне розв'язання. Невідому густину потенціалу простого шару 
обчислюють з нелінійного інтегрального фон баз тяжіння: 


1, с(х) сов(п, р) пе 2 
18 (х)-- У Ї ії с)48. ук в ее до соз( р, 4)45, 45. зах) хеду. () 








Иого розв'язок с), х є ду еквівалентний розв'язку задачі Алексідзе з граничними даними на 
поверхні ду Ляпунова, оскільки за будь-якого вибору густини потенціал простого шару задоволь- 
няє в області у" рівнянню Лапласа, а знайдене з (7) значення густини забезпечує виконання гра- 


ничної умови. Питання розв'язності задачі Алексідзе редукується до вияснення умов існування, 
єдиності 1 стійкості розв'язку рівняння (7). Це рівняння можна спростити до такого вигляду: 


)ч Пе І , | от). со5(р, 4)45, 45. -287(х), хеду: (8) 
гу ХО Х 


Розв'язки рівнянь (7, 3) допомагають визначати не лише потенціал У(х), хє у", ай значення 





МГИПСТ в будь-якій точці необмеженої області у з (1), або за зручнішою для обчислень форму- 
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лою ах) щ і; у Бас ро - | - ЯКа не вимагає обчислення двократного інтеграла. 


Зведення задачі Алексідзе з граничними даними на поверхні Ляпунова до розв'язання неліній- 
них інтегральних рівнянь сили тяжіння (7, 8) дозволяє на їх прикладі легко вивчити питання роз- 
в'язності її розв'язків, і ефективно знаходити чисельні наближені розв'язки для областей складної 
форми |Дубовенко, 2010). Для практичних обчислень послідовних наближень густини потенціалу 
простого шару за наведеним вище узагальненим алгоритмом доцільно застосувати такі ітераційні 
схеми: 


1. Со (х) ав «(х), О»0 (2) - 0, хе бу 7 Очі (2) 4 ор (х)ч ор (2) » о (2) З А, РИ со па 0,2 й. 


де А, Ід (х) -Боа,)- Гкибьеот дин (235, , К (обоє га а й а, 
бу Гі Хх- 


Мазедн ов") Всі) 


2. би о хєду, дині) тд о), он (х)2 Ас», (х) паз 0,0, 
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ду 
де суо(Х) з (2), сьо) 20, хєду. 


Розв'язання нелінійної граничної задачі Алексідзе для відновлення потенціалу за значеннями 
МГИПСТ означає відшукання потенціалу простого шару (4), густину якого відшукують з рівняння 
(5). Знаходити густину з еквівалентного нелінійного інтегрального рівняння сили тяжіння (7) не- 
доцільно, але воно зручне для вивчення умов коректності її постановки з граничними даними на 


класі поверхонь Ляпунова С ОУуг). На цьому класі задача Алексідзе редукована до розв'язання 
двох еквівалентних нелінійних інтегральних рівнянь (7, 8) сили тяжіння. У цих редукціях вона ко- 
ректна на парі банахових просторів В(у" 3) В(ду), до яких належать вхідні дані і шуканий розв'я- 


зок. 
Наступна стаття буде присвячена особливостям чисельних процедур за ітераційними схемами. 
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